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1 3.  ABSTRACT  (Maximum  200  words) 

A  study  is  underway  to  evaluate  the  proposed  ionic  mechanism  of  soot  formation  by  developing  a 
computer  model  and  comparing  it  with  experimental  results  and  with  equivalent  models  based  on 
neutral  species.  During  the  present  report  period  both  the  thermodynamic  and  reaction  kinetic 
data  base  for  ion-molecule  reactions  have  been  updated  and  expanded.  The  CHEMKJN-II  flame 
code,  revised  last  report  period,  to  handle  ions  is  now  operating  routinely  and  has  been  used  to 
calculate  ion  profiles  for  the  standard  acetylene/oxygen  flame.  The  results  are  generally  in  good 
agreement  with  experiment,  except  for  the  benzotropylium  ion,  for  which  the  calculated  value  is 
about  three  orders  of  magnitude  less  than  the  experimental  value.  Poor  thermodynamics  are 
suspected.  A  comparison  of  data  in  the  literature  on  the  effect  of  equivalence  ratio  on  acetylene 
and  soot  yield  with  values  calculated  by  thermodynamic  equilibrium  for  a  number  of  fuels  show 
that  in  general  measured  acetylene  and  soot  concentrations  are  far  in  excess  of  equilibrium, 
especially  below  the  threshold  soot  index.  The  standard  acetylene/oxygen  flame  seems  to  be  an 
exception.  The  results  continue  to  support  an  ionization  mechanism  of  soot  formation 
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I.  INTRODUCTION 


This  is  the  second  annual  report  on  a  program  to  develop  and  evaluate  a  model  for  an  ionic 
mechanism  of  soot  formation.  The  model’s  predictions  will  be  compared  with  experimental  data  to 
confirm  its  validity  or  demonstrate  its  inadequacy,  and  they  will  be  compared  with  corresponding 
predictions  of  models  based  on  neutral  mechanisms  of  soot  formation.  The  major  tasks  have  been  to 
develop  adequate  data  bases  of  ion  thermodynamics  and  ion-molecule  reaction  kinetics.  While 
thermodynamic  data  were  available  on  large  ions  with  even  numbers  of  carbon  atoms  there  were  no  such 
data  for  large,  odd  carbon  number  ions.  The  thermodynamic  data  base  proved  to  be  both  more  difficult 
to  prepare  and  also  more  important  than  the  ion  kinetic  data  base.  The  available  theory,  Langevin  theory, 
which  we  previously  used  to  calculate  ion-molecule  reactions  at  flame  temperatures,  was  inadequate  for 
large  ions.  We  extended  this  theory  to  cover  large  ions  and  found  that  a  slightly  positive  temperature 
coefficient  was  introduced  into  the  calculated  rate  constants.  However,  use  of  the  extended  theory  made 
very  little  difference  in  the  final  results.  The  thermodynamic  and  kinetic  data  bases  have  been  prepared 
and  are  continually  being  updated  as  new  data  appear  in  the  literature. 

In  the  last  report,1  we  described  revisions  to  the  Sandia  Flame  Code  to  enable  it  to  handle  ions. 
Problems  arose  with  the  modified  code  but  we  could  not  determine  whether  they  were  related  to  the  code 
itself  or  to  the  developed  mechanism.  Since  we  are  pioneering  a  new  area,  there  are  no  standards  against 
which  to  compare  our  results  for  confirmation  of  our  procedures.  During  the  last  year  we  have  increased 
our  confidence  in  the  computer  code  and  have  made  sufficient  improvements  in  the  thermodynamics  and 
kinetics  so  that  we  are  now  comfortable  with  the  results  and  feel  we  have  justified  the  effort  on  evaluating 
an  ionic  mechanism  of  soot  formation.  The  mechanism  holds  up  to  modeling  analysis  better  than  the 
neutral  mechanism,  despite  the  more  intensive  effort  that  has  gone  into  the  neutral  mechanism. 


We  initiated  an  effort  to  understand  the  relationship  of  soot  formation  and  the  thermodynamics 
of  soot  formation.  Experimental  data  in  the  literature  on  soot  threshold  equivalence  ratios,  soot  yields, 
and  acetylene  concentrations  are  compared  with  thermodynamic  equilibrium  calculations  of  these 
quantities.  The  results  are  surprising  and  vary  grossly  from  fuel  to  fuel. 
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II.  STATEMENT  OF  WORK 


The  ionic  mechanism  of  soot  formation  in  flames  will  be  further  evaluated  and  compared  with 
the  neutral  mechanism  by  pursuing  the  following  phases: 


Phase  I. 


1. 


2. 


3. 


4. 


Phase  II. 

1. 

2. 


3. 

Phase  III. 

1. 


Extend  the  Ionic  Mechanism  to  Benzene-Oxygen  Flames 

Extend  the  thermodynamic  and  reaction  rate  coefficient  data  base  to  include  those  species 
found  in  benzene/oxygen  flames  which  are  not  present  in  acetylene/oxygen  flames. 

Organize  the  experimental  data  available  on  the  benzene  flame  including  both  neutral  and 
ionic  species. 

Develop  a  detailed  ionic  mechanism  for  the  formation  of  soot  in  benzene  flames  and 
submit  this  to  others  to  run  on  a  large  computer  to  compare  its  agreement  with 
experimental  data  and  to  compare  its  simulation  of  soot  formation  with  the  neutral 
mechanism. 

Compare  the  computer  modeling  results  with  experimental  data  and  interpret  the  results 
in  terms  of  the  major  chemical  pathway  to  soot  in  the  benzene  flame. 

Model  Coagulation  and  Agglomeration 

Extend  the  detailed  mechanism  of  soot  nucleation  to  include  coagulation  of  large  ions  and 
neutrals,  and  charged  with  charged  and  neutral  incipient  particles. 

Program  a  desktop  computer  to  test  the  model  developed  in  Task  1  (in  a  limited  fashion) 
against  experimental  data.  If  warranted,  submit  this  model  to  others  to  test  on  a 
mainframe  computer. 

Interpret  the  results  from  Task  2  in  terms  of  a  simplified  mechanism. 

Develop  a  Theory  for  Large  Ion-Molecule  Reactions 

Extend  the  Langevin  theory  of  ion-molecule  reactions  to  include  large  ions  by  removing 
the  restriction  of  a  point  charge  on  the  ion. 
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Phase  IV. 


Compare  Thermodynamic  Predictions  of  Soot  Formation  with  Experimental  Observations 


1.  Collect  and  organize  the  literature  data  on  soot  yield  and  acetylene  concentrations  as  a 
function  of  equivalence  ratio. 

2.  Calculate  the  thermodynamic  equilibrium  concentration  of  acetylene  and  soot  as  a 
function  of  equivalence  ratio  for  the  literature  systems  identified  in  Task  1. 

3.  Compare  the  experimental  and  calculated  values  and  interpret  them  in  terms  of 
generalizations  relevant  to  the  mechanisms  of  soot  formation. 


III.  THERMODYNAMICS 

The  thermodynamics  of  cations  continues  to  be  a  major  source  of  uncertainty  in  the  kinetics 
model.  The  thermodynamics  of  several  cations  have  been  revised  according  to  the  most  recent  literature 
and  we  have  added  a  few  larger  cations  to  the  original  list  in  anticipation  of  incorporating  larger  ions  into 
the  kinetics  scheme.  As  before,  we  prefer  evaluated  literature  values  for  enthalpies  of  formation-; 
however,  when  they  are  not  available,  we  estimated  the  enthalpies  of  formation  from  proton  affinity, 
ionization  potential,  or  QM  calculated  values.  Heat  capacities  and  entropies  were  calculated  by  either 
rigid-rotor-harmonic  oscillator  statistical  mechanics  (SM)  formulae,3  hindered  rotor  harmonic  oscillator 
SM  formulae4  or  Benson’s  group  additivity  scheme5  using  the  most  recent  group  values.6  SM  formulae 
were  preferred  over  Benson’s  method  when  reliable,  experimental  or  high  level  QM,  structural  and 
vibrational  frequency  information  was  available  for  an  ion.  We  have  listed  the  changes  and  new  data,  in 
Table  I,  using  a  format  similar  to  that  of  Ref.  7. 

The  C^Hj4  (vinyl)  cation  is  dropped  from  our  original  species  list  because  recent  QM 
calculations  suggest  that  this  ion  is  a  transition  state  species  between  H-atom  shifted  cyclic 
structures.8  The  most  troublesome  species  in  our  kinetics  modeling  efforts  are  the  C,iH94  cations.  The 
most  recent  literature  values  for  enthalpies  of  formation  appear  to  be  consistent  with  our  earlier  estimates; 
heat  capacity  and  entropy  estimates  are  unchanged.  We  will  continue  to  search  for  more  stable  isomers 
of  this  ion. 

In  addition  to  the  cations,  we  have  estimated,  using  Benson  groups,  the  thermodynamics  of  a 
large  coronene  (8-circumcoronene)  with  the  molecular  formula  and  having  the  correct  C/H  ratio 

for  soot  (=  10).  We  utilized  this  species  in  our  equilibrium  calculations  for  determining  if  soot  is,  or 
is  not,  in  thermodynamic  equilibrium  in  flame  environments  (see  Section  VII). 
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IV.  REACTION  MECHANISM  AMD  REACTION  COEFFICIENTS 

During  this  period,  the  reaction  mechanism  has  remained  essentially  the  same,  with  some  changes 
in  isomers  used  and  the  deletion  of  others.  The  only  changes  in  rates  have  resulted  from  use  of  the 
extended  Langevin  theory  (results  to  be  submitted  for  publication  in  the  near  future)  to  calculate  the  rates 
of  ion-molecule  reactions. 

To  apply  the  extended  Langevin  theory  to  our  set  of  reactions,  we  empirically  reduced  it  to  two 
simple  equations  for  calculating  the  A  and  n  terms  in  the  classical  Arrhenius  equation,  (i.e.,  k(T)  =  A 
Tn  e-E3Ct/RT.  We  assumed  that  Eact  is  small  (—  0)  for  our  ion-molecule  reactions.  We  used  the 
conducting  sphere  model  of  the  extended  Langevin  theory.  For  acetylene  and  diacetylene,  the  two  major 
neutral  species  reacting  with  the  ions,  we  calculated  rate  coefficients  as  a  function  of  temperature  for  a 
series  of  ion  diameters.  The  polarizabilities  of  acetylene  and  diacetylene  (important  parameters  in  the 
calculations)  are:  3.33  x  10~"4  and  5.99  x  10~24  cm3,  respectively.  The  calculated  rate  coefficients 
for  each  ion  diameter  were  then  fit  by  regression  analysis  to  the  equation: 

log  k  =  log  A  +  n  log  T.  (1) 

This  gave  a  table  of  A’s  and  n’s  as  a  function  of  ion  diameter.  These  were  then  each  fit  to  simple 
equations  giving  the  dependence  of  A  and  n  on  ion  diameter  (d)  for  both  acetylene  and  diacetylene. 

Acetylene: 

A  =  (0.337  +  1.2  d)  x  1014,  cm3/mol/s 

n  =  0.292  +  0.0397  x  d,  d  in  nm 

Diacetylene: 

A  =  (2.96  +  9.93  x  d)  x  1013,  cm3/mol/s 

n  =  0.298  +  0.033  x  d,  d  in  nm 

For  all  of  the  fits,  the  correlation  coefficients  (r)  were  0.98  or  greater,  for  d  S  0.7  nm;  there  were  large 
changes  in  values  below  this  ion  diameter  which  is  where  the  extended  Langevin  theory  breaks  down. 

The  size  of  the  ions  was  deduced  from  experimental  ion  mobility  (*)  data  for  PAH  molecules  in 
nitrogen  gas.9  From  the  mobility  data  and  the  hard  sphere  diameter  of  the  nitrogen  molecule  (=  0.375 
nm.  Ref.  10),  using  Langevin’s  ion  mobility  equation,11 
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we  back  calculated  the  ion  diameters  of  the  PAH  ions  in  nitrogen  gas.  Here,  M  is  the  MW  of  nitrogen 
and  m  is  the  MW  of  the  ion,  K  is  the  dielectric  constant  of  nitrogen  and  p  is  the  gas  density.  A  is  a 
complicated  function  of  a  parameter  (X  -  not  mean  free  path),  gas  pressure,  p,  K  and  the  sum  of  the 
neutral  and  ion  collision  radii  (D12).  The  X  parameter  is  given  by: 

=  %*pP\2  (3) 

(*  -  1)  e2 

Fitting  these  data  to  various  empirical  functions,  we  found  the  following  relationship  between  ion  size 
and  ion  mass: 

ln(d.  A)  =  -1.375  +  2.154  In  (ln(m))  (4) 

The  correlation  coefficient  (r2)  for  this  fitting  equation  was  0.999-almost  a  perfect  fit  to  the  ion 
diameters. 

For  comparison,  we  plot  ion  diameters  deduced  from  ion  mobility  data  against  ion  diameters 
calculated  for  homogeneous  spheres  having  densities,  p,  of  1.5  or  2.0  g/cm3-the  density  range  for  soots. 
The  results  are  plotted  in  Fig.  1  and  they  show  close  agreement.  The  fitting  function,  Eq.  (4),  is 
routinely  used  to  calculate  PAH  ion  diameters  for  all  ions  in  our  modeling  work. 

The  mechanisms  for  ion-molecule  reactions  used  in  the  computer  runs  to  be  presented  here  are 
displayed  in  Table  II,  in  which  the  previous  reaction  coefficients  are  compared  with  those  calculated  by 
the  extended  Langevin  theory.  There  is  little  difference.  There  is  some  improvement  in  the  temperature 
coefficient  and  the  increased  rates  for  high  molecular  weight  species  improves  the  agreement  between 
calculated  and  experiment,  see  below. 

In  the  kinetic  scheme  presented  in  Table  II,  we  use  experimental  profiles  for  C3H3  +  and  H5C6+ 
(phenyl  cation)  because  we  do  not  have  good  experimental  profiles  for  the  HCO+  ion,  the  initial  ion 
produced  by  chemiionization,  from  which  all  other  ions  grow.  When  we  calculated  the  concentration  of 
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HCO  +  ,  using  as  our  basic  reaction  scheme  the  neutral  mechanism  for  soot  formation  and  adding  the 
chemiionization  reactions  to  form  HCO+,  we  calculated  a  value  too  high  compared  to  the  experimental 
value.  We  will  return  to  this  task  in  future  work. 

How  to  handle  C3H3  +  isomers,  propargylium  and  cyclopropenium  cations,  is  not  clear.  The  ion 
is  identified  as  a  flame  ion  by  mass  spectrometry  which,  unfortunately,  provides  no  direct  information 
on  the  ion  structure.  Two  isomeric  structures  for  C3H3+  have  been  considered  as  potentially  important. 
(These  are  the  linear  propargylium  cation,  denoted  as  C3H3(1)+,  and  the  cyclic  cyclopropenium  cation, 
C3H3(c)  +  .)  We  have  chosen  to  identify  the  C3H3(!)+  ion  with  the  experimentally  observed  ion  in  the 
mechanism.  Recent  calculations12  show  that  the  cyclic  cation  is  more  stable  than  the  linear  cation  by  1 16 
kj/mol;  the  experimental  difference  is  105  =  - 17  kj/mol.13  The  difference  in  stability  decreases  with 
increasing  temperature;  on  the  basis  of  the  derived  thermodynamics  at  300  K  the  equilibrium  ratio  of 
linear  to  cyclic  is  5  x  10~18  but  at  2000  K  the  ratio  is  0.03,  Fig.  2. 

It  has  been  determined  experimentally  that  at  low  temperatures  the  cyclic  isomer  is  not  reactive 
with  acetylene  or  diacetylene,  but  that  the  linear  isomer  is  reactive  with  acetylene  and  diacetylene. 14 
There  is  no  evidence  that  the  cyclopropenium  cation  isomerizes  to  the  propargylium  cation  in  the  absence 
of  encounter  complexes  with  neutral  molecules,  e  g.,  acetylene. 12  Cameron  et  al.12  suggest  that  in 
sooting  flames  the  cyclic  form  may  not  be  the  most  abundant.  In  the  flame,  the  C3H3+  cation  is  assumed 
to  be  formed  in  a  number  of  elementary  reactions  involving  HCO  +  .15  Based  on  the  rate  constants  and 
thermodynamics  we  adopted  in  developing  the  ionic  mechanism,  under  a  variety  of  experimental  flame 
conditions,  computer  experiments  showed  that  these  two  species  are  in  equilibrium.  This  problem 
requires  further  consideration. 


V.  COMPUTER  CODE 

The  revision  of  the  computer  code  to  accommodate  input  of  experimental  profiles  of  neutral 
species  involved  in  the  ionic  mechanism  has  made  it  much  easier  to  test  the  ionic  mechanism 
independently  of  the  neutral  species  model.  It  has  also  greatly  reduced  computer  time.  Considerable 
effort  was  committed  to  getting  the  modified  computer  code  to  work  consistently.  Some  of  the  difficulties 
we  reported  in  the  last  report1  were  due  to  imbalance  of  the  cation  and  electron  concentrations;  this  was 
corrected  in  the  program.  Some  of  the  other  problems  were  corrected  by  going  to  a  smaller  grid  size 
(i . e . ,  closer  spacing  of  calculated  concentrations  vs.  distance  steps).  The  program  is  now  used  routinely 
and  with  confidence. 
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VI.  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  ION 
CONCENT  RATION  PROFILES 


In  the  last  report1  we  noted  difficulties  we  were  having  getting  reasonable  agreement  between 
experimental  and  calculated  ion  profiles.  In  fact,  since  the  inception  of  this  program  there  has  been  a 
major  problem  with  our  model  in  simulating  the  shape  and  position  of  the  ion  profiles  with  respect  to  the  I 

experimental  profiles.  We  were  also  experiencing  problems  with  modeling  the  observed  double  peaked 
ion  profiles.  There  are  frequently  two  peaks  in  the  experimental  ion  profiles,  the  peak  closest  to  the 
burner  being  relatively  small.  The  first  peak  in  the  model  appeared  to  be  due  to  the  ion-molecule 
reactions  having  large  entropy  terms  so  that  at  the  lower  flame  temperatures  associated  with  the  position 
of  the  first  peak  the  forward  rates  dominated.  We  considered  many  possible  causes.  We  found,  for  * 

example,  that  assigning  activation  energies  to  ion-molecule  reactions  helped.  We  even  went  so  far  as  to 
consider  abandoning  the  hypothesis  that  ions  play  an  important  role  in  soot  formation,  seeking  other 
explanations  of  the  experimental  data. 

We  are  thus  pleased  to  report  that  without  any  drastic  treatment  of  the  data  or  without  making 
any  ad  hoc  assumptions  we  have  obtained  a  very  satisfactory  agreement  between  the  computer  model  and 
experiment.  Computed  and  observed  ion  profiles  are  compared  in  Figs.  3-5  for  the  standard 
acetylene/oxygen  flame  at  2.67  kPa  and  equivalence  ratio  =  3.0.  The  calculated  ion  profiles  were  made 
using  the  kinetic  scheme  in  Table  II.  with  the  rate  coefficients  calculated  by  the  extended  Langevin  theory 
for  conducting  spheres  and  reverse  reactions  calculated  with  the  latest  thermodynamics.  While  a  number 
of  isomers  are  included  in  the  calculations,  only  the  most  abundant  cation  isomer  for  any  ion  mass  is 
plotted;  the  mass  spectrometer  is  blind  to  isomeric  structure.  It  is,  of  course,  conceivable  that  the  most 
thermodynamically  stable  isomer  cannot  be  formed  because  of  reaction  kinetic  restrictions.  This  will  have  I 

to  be  pursued  subsequently;  if  we  can  simulate  the  experiment  with  reasonable  assumptions,  we  will 
examine  the  assumptions  in  more  detail. 

In  the  computer  model:  the  experimental  cation  profiles  for  C3H3  +  (assumed  to  be  propargyl, 
see  above),  and  H5C6+  (stable  cyclic  isomer)  were  used  as  input.  These  are  the  lowest  mass  ions  in  this 
flame  for  which  the  experimental  data  are  considered  reliable  or  are  available.  Because  the  ion  growth 
paths  seem  to  run  parallel  for  even  and  odd  numbers  of  carbon  atoms,  and  generally  independent  of  one 
another,  it  was  necessary  to  use  both  an  odd  and  an  even  numbered  carbon  cation  as  input.  However, 
there  are  a  few  ion-molecule  reactions  with  C3H4  in  our  model,  and  these  reactions  couple  the  odd  and 
even  ion  growth  pathways.  All  neutral  species  profiles  that  enter  into  the  mechanism  (i.e.,  C2H2,  C4H2. 

C3H4  and  H2)  were  also  used  as  input  All  of  the  experimentally  observed  ions  are  included  in  the 
mechanism  and  no  species  are  included  that  have  not  been  observed.  This  is  an  important  point  when 
comparing  these  results  with  results  of  the  neutral  mechanism  models,  because  in  application  of  those 


7 


TP-518 


models  many  unobserved  species  are  included  and  some  of  the  observed  neutral  species  are  not 
calculated.  No  changes  have  been  made  in  the  ion-electron  recombination  reactions  or  their  rates  since 
the  last  report.  Negative  ion  formation  and  recombination  have  not  been  included.  Experimental  profiles 
were  taken  from  the  literature. 

A  temperature  profile  was  also  used  as  input;  it  was  generated  from  the  experimental  profile  by 
keeping  the  temperature  at  the  burner  fixed  at  the  experimental  temperature  and  scaling  the  rest  of  the 
experimental  profile  by  the  formula: 


1  scaled 


=  T. 


exp 


1-0.1 


exp 


990 


2020  -  990 


(5) 


where  Texp  =  original  experimental  value 

2020  =  original  maximum  experimental  flame  temperature 

990  =  extrapolated  experimental  flame  temperature  at  0.91  mm  above  burner  surface 

This  reduced  the  maximum  flame  temperature  to  1818  K.  This  is  within  the  accuracy  of  the 
thermocouple  temperature  measurements  and  in  the  expected  direction. 

From  Figs.  3-5  it  is  obvious  that  the  shapes  and,  in  general,  the  concentration  levels  of  the 
experimental  profiles  are  well  simulated  by  the  model.  Clearly  there  still  remains  a  question  with  respect 
to  the  CnH9  +  ion  which  carries  over  to  the  two  largest  cations,  Ci;;H9+  and  C13H9+,  Fig.  5.  Of 
concern  is  the  fact  that  some  of  the  calculated  concentration  profiles  show  an  increase  as  the  burner  is 
approached.  As  noted  above,  double  peaks  have  been  observed  experimentally6  but  the  effect  here  is 
greater  than  in  the  experiments  and  there  seems  to  be  no  particular  correlation  with  ion  size,  odd/even 
carbocations  or  the  thermodynamic  stability  of  the  ion  In  addition,  the  experimental  results  are  very 
limited  close  to  the  burner  where  there  are  questions  about  interference  of  the  mass  spectrometer  sampling 
cone  with  the  flame.  We  feel  that  the  calculated  increase  toward  the  burner  is  more  an  artifact  of  the 
model  which  has  too  great  a  temperature  effect  (thermodynamics)  and  uses  flame  temperature  profiles 
that  are  not  that  accurate  near  the  burner  port.  This  situation  could  be  corrected  by  incorporating 
activation  energies  into  some  of  the  ion-molecule  reactions  and/or  modifying  the  input  temperature  profile 
even  more  than  we  have  already  done. 

To  better  evaluate  the  agreement  between  experimental  and  calculated  concentrations,  the 
maximum  concentration  ratios,  log l0  (experimental  concentration/calculated  concentration)  are  presented 
in  Fig  6.  The  isomer  in  closest  agreement  with  experiment  is  plotted  and  the  calculated  results  using 
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Langevin  and  extended  Langevin  are  compared.  A  bar  to  the  right  of  zero  indicates  that  the  experimental 
concentration  exceeds  the  calculated  value  and  to  the  left,  the  opposite.  The  source  of  disagreement  with 
the  CnH9  +  cation  has  been  considered  by  carefully  reviewing  the  thermodynamics,  including  the 
possibility  of  other  more  stable  isomers,  without  arriving  at  a  solution.  It  appears  that  there  is  little 
latitude  for  adjusting  the  rate  coefficients.  The  experimental  results  could  be  called  into  question,  but  we 
feel  that  the  major  error  remains  in  the  derived  thermodynamics.  This  will  be  pursued  further  as  we 
extend  the  model  to  larger  ions;  development  of  the  thermodynamics  and  reaction  kinetics  for  these  ions 
has  already  been  completed. 

Oxidation  of  ions  has  yet  to  be  added  to  the  mechanism.  Unlike  the  neutral  mechanism,  wherein 
oxidation  removes  the  neutral  growth  species  from  the  growth  chain,  for  ions,  oxidation  simply  removes 
carbon  atoms  from  the  ion,  making  a  smaller  ion,  or  removes  hydrogen  from  the  growing  ion  at  a  greater 
rate.  Removal  of  hydrogen  may  in  fact  promote  ion  growth  to  soot  which  has  a  C/H  ratio  of  about  10/1, 
while  for  the  largest  ion  treated  here  the  ratio  is  only  1.4/1.  Since  the  ions  .onsidered  in  the  model  have 
a  decreasing  free  energy  per  carbon  atom  with  increasing  size,39  it  is  assumed  that  the  charge  will  remain 
on  the  larger  hydrocarbon  products,  thus  e.g . : 

C13H9+  +  OH  C13H8+  +  H20  (6) 

would  he  expected  to  occur,  rather  than: 

C13H9  +  +  OH  -*  C13H7  +  H30+  (7) 

This  of  course  requires  further  consideration  with  other  oxidizers  and  other  products. 

VII,  THERMODYNAMICS  OF  SOOT  FORMATION 

The  starting  point  in  understanding  any  chemical  system  is  the  thermodynamics  of  the  system. 
This  has  not  been  done  for  the  soot  forming  problem  to  any  great  extent.  Millikan17  was  probably  the 
first  to  point  out  that  the  critical  equivalence  ratio  at  which  soot  first  appears  is  not  consistent  with 
equilibrium  theory.  This  point  has  been  repeated18-37  but  in  all  cases  the  thermodynamics  was  limited 
to  assuming  soot  was  equivalent  to  graphite.  Soot  in  fact  is  about  10  mol  percent  hydrogen  and  there  are 
no  thermodynamic  data  on  real  soot.  We  have  thus  initiated  a  study  to  evaluate  the  thermodynamics  of 
the  soot  formation  process  which  includes  estimating  the  thermodynamics  of  a  typical  soot  and  including 
in  the  equilibrium  calculation  large  polycyclic  aromatic  hydrocarbons  which  have  previously  been 
neglected  in  such  calculations.  In  a  previous  publication  we  showed  that  the  thermodynamic  driving  force 
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to  produce  large  hydrocarbon  species  was  greater  for  ions  than  for  free  radicals  up  to  a  mass  of  about 
350  —400  u  where  they  become  equal.1 

The  motivation  for  our  interest  is  twofold.  In  some  earlier  work  at  AeroChem19'20  it  was 
observed  that  the  threshold  soot  index,  TS1,  for  a  number  of  fuels  occurred  at  a  fixed  experimental 
temperature  which  was  much  lower  than  the  adiabatic  flame  temperature.  The  implication  is  that  the 
energy  which  should  have  gone  into  raising  the  gas  temperature  was  expended  in  making  soot.  It  seems 
unrealistic  to  assume  that  radiation  was  the  major  factor  since  the  soot  loadings  are  small  at  threshold. 
Another  motivation  was  to  better  understand  a  problem  associated  with  the  ionic  mechanism  of  soot 
formation.  In  a  fuel  rich  system  the  total  ion  concentration  decreases  with  increasing  equivalence  ratio 
until  at  some  critical  equivalence  ratio  beyond  the  soot  threshold  the  ion  concentration  again  increases 
and  negative  ions  are  observed.21  At  the  threshold  the  ion  composition  changes  and  large  ions  are 
observed.  The  question  is,  why  does  soot  appear  when  the  ion  concentration  is  falling  if  ions  are 
important  to  the  mechanism  of  soot  formation?  Are  the  large  ions  formed  from  soot  or  is  soot  formed 
from  the  large  ions,  or  is  the  simultaneous  appearance  of  soot  and  large  ions  an  accidental  occurrence? 
Large  ion  formation  from  soot  particles  has  been  eliminated.22 

In  approaching  this  problem  we  collected  literature  data  on  soot  yield  and  acetylene  concentrations 
and  measured  flame  temperatures  as  a  function  of  equivalence  ratio.  We  then  calculated  these  quantities 
using  an  equilibrium  program.  We  added  to  the  thermal  data  base  of  that  program,  the  thermodynamics 
of  large  aromatic  hydrocarbons,  see  Table  III.  For  the  initial  calculations  we  assumed  the  soot  was 
graphite.  The  results  for  a  number  of  fuel  systems  are  shown  in  Figs.  7-13.  The  equivalence  ratio  at 
which  the  product  of  a  balanced  reaction  would  be  CO  +  H2  is  indicated  on  the  figure.  Except  for 
acetylene,  there  is  a  dramatic  difference  between  the  experimental  and  calculated  concentrations  of 
acetylene  and  soot  and  where  they  appear  with  respect  to  equivalence  ratio.  They  clearly  appear  at 
equivalence  ratios  much  lower  than  where  they  are  predicted  to  occur  thermodynamically.  The 
differences  in  the  concentrations  of  calculated  and  observed  acetylene  are  especially  dramatic,  varying 
from  a  factor  of  10  for  acetylene  to  108  for  methane.  It  is  unfortunate  that  the  standard  low  pressure 
acetylene/oxygen  flame  which  has  received  so  much  attention  seems  to  behave  so  differently  from  other 
hydrocarbons  (Fig.  12). 


To  obtain  an  estimate  of  the  thermodynamics  of  soot,  we  examined  the  C/H  ratios  for  the  fully 
condensed  benzene  ring  compounds  of  the  circumcoronene  series  (C6N.NH6N).  At  N  =  10,  the 
molecular  formula  is  C600H60  and  this  molecule  has  the  same  C/H  ratio  as  soot  (~  10).  Using  Benson's 
group  additivity  method5,6  to  estimate  the  thermodynamic  quantities  (i.e  ,  AH°(298  K),  S°(298  K)  and 
Cp(T)),  we  derived  the  thermodynamics  of  this  gas  phase  molecule  (Table  I).  To  reduce  the  gas  phase 
data  appropriate  to  solid  phase  "soot",  we  required  the  enthalpy  and  entropy  of  sublimation,  and  to  get 
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these  values,  we  had  to  estimate  C600H60's  boiling  point  (bp).  Using  experimental  boiling  points  for 
benzene,  naphthalene,  phenamhrene,  and  coronene,  we  modified  Joback’s  estimation  method34  to  give 
close  agreement  with  the  experimental  boiling  points.  Then  we  used  this  method  to  obtain  the  bp  of 
CeooHgo-  From  the  estimated  bp,  Vetere's  relationship35  was  used  to  obtain  the  enthalpy  and  entropy 
of  sublimation  for  These  values  were  then  subtracted  from  the  gas  phase  values  to  give  the 

enthalpy  of  formation  and  entropy  of  solid  C^ooH^. 

When  this  derived  thermodynamic  data  tor  "soot”  was  incorporated  into  the  adiabatic  equilibrium 
calculation,  there  was  no  significant  difference  in  calculated  soot  yields  nor  the  equivalence  ratio  at  which 
soot  appears.  The  threshold  equivalence  ratio  for  soot  formation  seems  to  depend  more  upon 
stoichiometry  than  thermodynamics.  Note  in  Figs.  7  through  13  the  indication  on  the  equivalence  ratio 
axis  of  where  the  products  are  determined  by  stoichiometry  to  be  CO  +  H2.  Above  this  composition, 
increases  in  the  fuel  should  lead  to  increased  yields  of  hydrocarbons,  due  to  insufficient  oxygen. 

Since  soot  is  a  supra  equilibrium  product  it  requires  a  driving  force  from  a  nonequilibrium 
precursor.  Something  must  be  present  in  excess  of  equilibrium  to  initiate  the  process.  For  the  free 
radical  mechanism,  it  is  assumed  that  the  hydrogen  atom  concentration  is  in  excess  of  equilibrium  and 
is  the  driving  force.36  For  the  ionic  mechanism  it  is  assumed  that  chemiions,  formed  by  a  nonequilibrium 
process  and  roughly  two  orders  of  magnitude  in  excess  of  equilibrium,37  are  the  driving  force.  The 
importance  of  nonequilibrium  acetylene,  and  presumably  diacetylene,  concentrations  seems  apparent  from 
inspection  of  Figs.  7  to  13.  The  hydrogen  atom  or  ion  might  still  play  a  role,  because  it  has  been 
demonstrated  that  polymerization  of  acetylene  alone  is  not  the  mechanism38;  otherwise  polyacetylenes 
would  be  the  precursor  species  and  it  is  generally  agreed  that  the  precursor  species  is  a  polycyclic 
compound,  ion,  or  neutral.  Since  more  acetylene  is  involved  in  the  growth  phase  of  the  ionic  mechanism 
than  ions,  only  one  ion  per  carbon  particle  is  involved  and  the  global  rate  of  soot  formation  may  be 
written: 


=  k  x  /  x  An  dt  (8) 

dt 


where  k  is  a  rate  constant,  I  is  the  concentration  of  ions,  A  is  the  concentration  of  acetylene  and  n  is  an 
integer  very  much  greater  than  1 .  Thus,  both  ions  and  acetylene  must  be  present,  and  this  occurs  at  an 
equivalence  ratio  where  the  ion  concentration  is  falling  and  the  acetylene  concentration  is  rising  with 
increasing  equivalence  ratio. 
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VIII,  PRESENTATION 


Calcote,  H.F.,  "Modeling  Study  to  Evaluate  the  Ionic  Mechanism  of  Soot  Formation,"  presented  at 
AFOSR  Contractors  Meeting  on  Propulsion  Research,  La  Jolla,  CA,  15-19  June  1992. 
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collecting  and  calculating  the  thermodynamics  data  and  for  running  the  modified  Sandia  Flame  Code. 
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TABLE  I 

HEAT  OF  FORMATION  AT  298  K,  ENTROPY  AT  298  K  AND  HEAT  CAPACITY  FOR 
CATIONS  AND  8-CIRCUMCORONENE  AT  1  ATMOSPHERE  PRESSURE  IN  THE  GAS  PHASE 

AfH°  units  are  kj/mol.  S°  and  Cp°  units  are  J/moI/K.  The  first  column  lists  the  molecular  formula 
below  which  is  the  name  used  in  the  kinetics  model.  GA  in  "Ref."  column  indicates  thermodynamic 
quantities  were  all  estimated  using  Benson’s  group  additivity  method. 

C„  AT  TEMPERATURE,  K 


Name 

AfH! 

s° 

298 

m. 

1000 

1500 

2000 

2500 

CoH  3  + 
C2H3  + 

1100 

221 

46.5 

60.4 

80.7 

91.6 

97.5 

100.9 

C4H3  + 
C4H3  + 

1217 

269 

70.2 

92.9 

123.6 

138.2 

145.5 

149.6 

C4H5+(a) 
C4H5  + 

1029 

290 

77.6 

111.5 

154.7 

175.1 

185.6 

191.5 

C4H5+(b) 

H5C4+ 

992 

278 

73.2 

109.1 

156.6 

177.1 

187.3 

192.9 

C5H3+(a) 
C5H2H  + 

1376 

295 

84.1 

110.1 

144.8 

161.0 

169.2 

173.7 

C5H3+(b) 
H3C5  + 

1301 

288 

78.6 

108.2 

144.6 

160.8 

169.0 

173.5 

C6H5+(a) 
H5C6  + 

1127 

282 

79.1 

129.9 

194.7 

222.1 

235.5 

242.4 

C6H5+(b) 
C6H5  + 

1320 

338 

107.5 

150.1 

199.4 

222.0 

233.8 

240.0 

C7H5  + 
C7H4H  + 

1237 

340 

90.0 

145.1 

214.2 

242.9 

256.9 

264.1 

C,,H9+(a) 
C11H9  + 

1054 

385 

157.7 

249.0 

365.1 

414.8 

438.4 

451.7 

16 

ISjI 
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Cp  AT  TEMPERATURE,  K 


Name 

AfHI 

S° 

298 

JQQ_ 

1000 

1500 

2000 

2500 

Ref. 

C„H9+(b) 

H9C10+ 

1051 

403 

150.9 

232.5 

339.6 

391.2 

422.5 

443.0 

12 

CuH,+(c) 
C11H8H  + 

962 

386 

148.5 

243.1 

363.4 

414.7 

438.8 

451.5 

3 

CI3H9+ 

C13H6H3+ 

918 

364 

169.4 

276.9 

409.5 

463.7 

489.7 

503.0 

13 

C54H19  + 

C54H19  + 

1314 

859 

613.7 

1041.7 

1476.3 

1621.5 

1687.4 

1725.8 

14 

c96h25+ 
C96H25  + 

1701 

1251 

1020.3 

1774.0 

2501.0 

2729.1 

2840.8 

2926.2 

14 

4880 

5436 

5532 

10244 

14276 

15316 

15759 

16026 

GA 

C600H60 
"gas  phase 
soot" 
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TABLE  II 

REACTIONS  AND  RATE  COEFFICIENTS  (LANGEVIN  AND  EXTENDED  LANGEVIN)  FOR 
IONIC  MECHANISM  OF  SOOT  FORMATION 


k(T)  =  AT"  exp  (Elct/RT) 
units  J,  cm3  mole,  s 


REACTIONS 

Langevin  Rate 

Parameters 

A  n  Eact 

Extended  Langevin  Rate 

Parameters 

A  n  Eact 

c3h3  + 

+ 

c2h2 

c5h2h+ 

4 

H, 

6.50e+  14 

0 

0 

8.60e+ 13 

0.3 

0 

c3h3+ 

+ 

c2h2 

- 

C5Hs  + 

6.50e+  14 

0 

0 

8.60e+  13 

0.3 

0 

c3h3" 

+ 

c4h2 

= 

c5h2h+ 

4 

CiH2 

7.40e+  14 

0 

0 

7.28e+  13 

0.3 

0 

c3h3+ 

+ 

c4h2 

= 

c7h4h+ 

7.40e+  14 

0 

0 

7.28e+  13 

0.3 

0 

c4h5+ 

+ 

c2h2 

= 

h5c6+ 

4 

h2 

6.20e+  14 

0 

0 

9. 17e+  13 

0.3 

0 

CsH2H  + 

+ 

H2 

= 

c5h5+ 

9.00e+  14 

0 

0 

8.97e+  14 

0.0 

0 

c5h2h+ 

4* 

c2h2 

= 

c7h4h+ 

6.00e+  14 

0 

0 

9.51e+  13 

0.3 

0 

CsH2H+ 

+ 

c3h4 

= 

h5c6+ 

4 

C,H, 

6.50e+  14 

0 

0 

9.51e+  13 

0.3 

0 

c5h5+ 

+ 

c2h2 

= 

c7h4h+ 

4 

H, 

6.00e+  14 

0 

0 

9.58e+  13 

0.3 

0 

c5h5+ 

+ 

c2h2 

= 

c7h5+ 

4 

H2 

2.40e+  14 

0 

-319 

9.58e+  13 

0.3 

0 

c5h5+ 

+ 

c4h2 

= 

c9h7+ 

6.50e+  14 

0 

0 

9.58e+  13 

0.3 

0 

c5h5+ 

+ 

c4h2 

= 

h7c9+ 

6.50e  +  14 

0 

0 

8.09e+  13 

0.3 

0 

c5h5+ 

+ 

c2h2 

c7h7+ 

6.50e+  14 

0 

0 

8.09e+  13 

0.3 

0 

h5c6+ 

+ 

c2h2 

= 

csh7+ 

5.80e+  14 

0 

0 

9.94e+  13 

0.3 

0 

c7h4h+ 

+ 

H2 

= 

c7h7+ 

8.90e+  14 

0 

0 

8.93e+  14 

0.0 

0 

c7h4h+ 

+ 

h2 

= 

h7c7+ 

8 . 90e  +  1 4 

0 

0 

8.93e  +  14 

0.0 

0 

c7h5  f 

4 

c2h2 

= 

c9h7+ 

5.70e+  14 

0 

0 

I.03e+ 14 

0.3 

0 

c7h5+ 

4 

c2h2 

= 

h7c9+ 

5.70e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h5+ 

4 

M 

= 

c7h4h+ 

4 

M 

4.60e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h5+ 

4 

c3h4 

= 

c8h7" 

4 

CiH2 

6. 15e  +  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h5  + 

4 

c3h4 

= 

c10h9+ 

6. 15e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h4h+ 

4 

c2h2 

h7c9+ 

5.70e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h4h+ 

4 

c2h2 

= 

c9h7+ 

5.70e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h7+ 

4 

M 

= 

H  7C7  + 

4 

M 

4.60e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h7+ 

4 

c2h2 

= 

h7c/ 

4 

H2 

5.70e+  14 

0 

0 

1.03e+  14 

0.3 

0 

h7c7+ 

4 

c2h2 

= 

h7c9+ 

4 

h2 

5.70e+  14 

0 

0 

1.03e+  14 

0.3 

0 

c7h7+ 

4 

c4h2 

= 

CuH8H  + 

6.10e+  14 

0 

0 

8,69e  +  13 

0.3 

0 

h7c/ 

4 

c4h2 

= 

C„HgH  + 

6. 10e+  14 

0 

0 

8.69e+ 13 

0.3 

0 
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Langevin  Rate  Extended  Langevin  Rate 


Parameters  Parameters 


REACTIONS 


c„h7+ 

+ 

c2h2 

= 

c10h9+ 

c8h7+ 

+ 

c2h2 

= 

H9C,o+ 

c8h7+ 

+ 

c3h4 

= 

ChH8H+ 

+ 

H2 

c8h7+ 

+ 

c4h2 

= 

c12h9+ 

c9h7+ 

+ 

M 

= 

h7c9+ 

+ 

M 

c9h7+ 

+ 

c2h2 

= 

CuH8H  + 

h7c9+ 

+ 

c2h2 

= 

c„h8h+ 

c9h7+ 

+ 

c4h2 

= 

ci3h6h3+ 

h7c9+ 

+ 

c4h2 

= 

ci3h6h3+ 

h9c10* 

+ 

M 

c10h9+ 

+ 

M 

C10H,+ 

+ 

c2h2 

= 

c12h9+ 

H2 

H0C10  + 

+ 

c2h2 

= 

c12h9+ 

h2 

h9c10+ 

+ 

c2h2 

= 

h9c12+ 

+ 

H2 

c10h9+ 

c3h4 

= 

ci3h6h3+ 

+ 

2H2 

c10h9  + 

+ 

c4h2 

= 

c12h9+ 

+ 

c2h2 

H9C,o  + 

+ 

C4H2 

= 

c12h9+ 

+ 

c2h2 

H9C,o* 

+ 

c4h2 

= 

h9c12+ 

+ 

C2H2 

c,,h8h  + 

+ 

C2H2 

= 

ci3h6h3+ 

+ 

H*2 

Ci,H9* 

+ 

c,h2 

= 

C13H6H3  + 

+ 

H2 

H,c,r 

+ 

M 

= 

c,,h8h+ 

+ 

M 

h9c„  + 

+ 

M 

= 

cuh9+ 

+ 

M 

c„h9  + 

+ 

M 

c,,h8h+ 

+ 

M 

h9c12+ 

+ 

M 

= 

c,;h9+ 

+ 

M 

c3h3  + 

+ 

E 

- 

c2n2 

+ 

CH 

c4h5  + 

+ 

E 

— * 

c,h2 

+ 

c2h3 

c5h2h+ 

+ 

E 

- 

c3h3 

+ 

C2 

c5h5  + 

+ 

E 

— * 

c3h3 

+ 

h5c6+ 

+ 

E 

-* 

c4h4 

+ 

c2h 

c7h4h  + 

+ 

E 

-* 

c4h2 

+ 

C3H3 

h7c7+ 

+ 

E 

-* 

c6h4 

+ 

ch3 

c8h7+ 

+ 

E 

—* 

c«h6 

+ 

c2h 

»¥V 

+ 

E 

— * 

c8h6 

+ 

CH 

c,0h9+ 

+ 

E 

CioH, 

+ 

H 

c„h8h+ 

+ 

E 

— * 

C,oH* 

+ 

CH 

A 

n 

Eact 

A 

n 

Eact 

5.64e+  14 

0 

0 

1.06e+  14 

0.3 

0 

5.64e+ 14 

0 

0 

1.06e  + 14 

0.3 

0 

6.02e+  14 

0 

0 

1.06e+  14 

0.3 

0 

5.95e+  14 

0 

0 

8.94e+ 13 

0.3 

0 

4.50e+  14 

0 

0 

1.09e+  14 

0.3 

0 

5.60e+  14 

0 

0 

1.09e+  14 

0.3 

0 

5.60e+  14 

0 

0 

1.09e+  14 

0.3 

0 

5.80e+  14 

0 

0 

9. 16e+  13 

0.3 

0 

5.80e+  14 

0 

0 

9.16e+  13 

0.3 

0 

4.50e+  14 

0 

0 

1.12e+ 14 

0.3 

0 

5.53e+  14 

0 

0 

1. 12e+  14 

0.3 

0 

5.53e+  14 

0 

0 

1.12e+  14 

0.3 

0 

5.53e+  14 

0 

0 

1.12e+  14 

0.3 

0 

5.85e  + 14 

0 

0 

1. 12e+  14 

0.3 

0 

5.75e+  14 

0 

0 

9.40e+  13 

0.3 

0 

5.75e+  14 

0 

0 

9.40e+  13 

0.3 

0 

5.75e+  14 

0 

0 

9.40e+  13 

0.3 

0 

5.50e+  14 

0 

0 

1.14e+  14 

0.3 

0 

5.50e+  14 

0 

0 

1.14e+I4 

0.3 

0 

4.50e+  14 

0 

0 

1.14e+14 

0.3 

0 

4.50e+  14 

0 

0 

1 .14e  + 14 

0.3 

0 

4.50e  + 14 

0 

0 

1 . 1 4e  +  1 4 

0.3 

0 

4.47e t 14 

0 

0 

1.16e+  14 

0.3 

0 

1 .05e+  19 

-0.5 

0 

1.05e+  19 

-0.5 

0 

1 . 20e  +  1 9 

-0.5 

0 

1.20e+  19 

-0.5 

0 

I.33e+  19 

-0.5 

0 

1.33e+  19 

-0.5 

0 

1.40e+  19 

-0.5 

0 

1.40e  +  19 

-0.5 

0 

I.42e  +  19 

-0.5 

0 

1.42e  +  19 

-0.5 

0 

!  50e  + 19 

-0.5 

0 

1  50e  + 19 

-0.5 

0 

1.52e+  19 

-0.5 

0 

1 .52e  +  19 

-0.5 

0 

1.61e  +  19 

-0.5 

0 

1.61e+  19 

-0.5 

0 

1.70e+  19 

-0.5 

0 

1 ,70e+  19 

-0.5 

0 

1.79e+  19 

-0.5 

0 

1.79e+  19 

-0.5 

0 

I.86e+  19 

-0.5 

0 

1.86e  ¥  19 

-0.5 

0 
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Langevin  Rate 
Parameters 

_ RF.  ACTIONS _  A  n  Eact 

C,,H9+  +  E  -  CnHg  +  H  1 ,93e+ 19  -0.5  0 

c13h6h3+  +  E  -  c13h9  1.99e+ 19  -0.5  0 


Extended  Langevin  Rate 
Parameters 
A  n  Eact 

1 ,93e  + 19  -0.5  0 

1.99e+ 19  -0.5  0 
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TABLE  III 

SPECIES  ADDED  TO  THERMODYNAMIC  DATABASE 
OF  EQUILIBRIUM  PROGRAM 


diacetylene 

C4H2 

triacetylene 

c6h2 

benzene 

c6h6 

toluene 

c7h8 

quatracetylene 

c8h. 

azulcne 

c,0h8 

naphthalene 

C,oH8 

acenaphthylene 

c,2hs 

acenaphthane 

Cr.H.o 

phenalene 

C13H10 

phenanthrene 

C|4H10 

acephenanthrylene 

C16H10 

cyclopenta[cd)pyrene 

C18H10 

benzo[cd]cyclopenta[mn]pyrene 

1 1 

benzo[ghi]perylene 

CnH  p 

indeno|5,6,7, 1  -pqra|perylene 

C;4Hi; 

coronene 

C24H,2 

benzo[a]coronene 

0:8^14 

naphtho[8, 1 ,2-abc]coronene 

^3(^14 

ovalene 

C32H,4 

circumcoronene 

n 

t 

X 

oo 

buckminsterfullerene 

^60 

DIAMETER, 


FIGURE  I  PAH  ION  DIAMETERS 

PAH  ions  diameters  calculated  from  experimental  ion  mobility  (Langevin), 
and  spheres  of  density  1.5  cm3  (p=1.5)  and  2.0  cm3  (p  =  2.0). 


93-7 


FIGURE  2  LOG,0  OF  EQUILIBRIUM  CONSTANT  FOR  SYSTEM 
CjHj+(1)  =  C^c) 


DISTANCE,  cm 


FIGURE  5  CATION  CONCENTRATION  PROFILES 

Acetylene/Oxygen  Flames  at  2.67  kPa,  equivalence  ratio  =  3.0. 

Experimental  values  represented  by  broken  line;  calculated  values  represented  by  solid  line. 
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- ! - 1  I 

C5H2H  + 

— 

C7H4H  + 

- 

& 

C7H7  + 

- 

f 

C8H7  + 

_ g 

- 

H7C9  + 

I 

i 

C10H9  + 

777777] 

CT1H8H  + 

- 

V///////S  '/'////// A 

C12H9  + 

- 

'wUSS^ 

C13H6H3+ 

- 

ill 

- 

1  i 

LOGiq 

D  1  2  3  4 

(EXPTL.  CONC./CALC.  CONC.) 

FIGURE  6  COMPARISON  OF  MAXIMUM  CALCULATED  ION  CONCENTRATIONS 
USING  LANGEVIN  RATE  COEFFICIENTS  (SHADED)  AND  EXTENDED  LANGEVIN 
RATE  COEFFICIENTS  (SOLID)  TO  THE  MAXIMUM  EXPERIMENTAL  ION 

CONCENTRATIONS 

Acetylene/Oxygen  Flames  at  2.67  kPa,  equivalence  ratio  =  3.0. 
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FIGURE  7  EXPERIMENT  VS.  EQUILIBRIUM 
METHANE/OXYGEN  FLAMES,  101  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results:  A  -  C^Hj,  Ref.  23;  ♦  -  Soot,  AeroChcm; 

•  -  Soot,  Ref.  24;  A  -  Soot,  Ref.  25;  O-  Temperature,  Refs.  23,  24. 


EQUIVALENCE  RATIO 

FIGURE  8  EXPERIMENT  VS.  EQUILIBRIUM 
PROPANE/OXYGEN  FLAMES,  15  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results  (Ref.  26):  A  -  CjH^,  ■  -  Soot,  O-  Temperature 
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FIGURE  10  EXPERIMENT  VS.  EQUILIBRIUM 
BF.NZENE/ARGON/OXYGEN  FLAMES,  5.3  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results  (Ref.  29):  A  -  C^H?;  ■  -  Soot;  O-  Temperature. 


I 


FIGURE  11  EXPERIMENT  VS.  EQUILIBRIUM 
TOLUENE/AIR  FLAMES,  101  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results  (Ref.  20):  •  -  Soot,  N2-rich  air;  A  -  Soot,  air;  ■  -  Soot,  02-rich  air; 
O-  Temperature,  N2-rich  air;  A  -  Temperature,  air;  □  -  Temperature,  Q2-rich  air. 
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FIGURE  12  EXPERIMENT  VS  EQUILIBRIUM 
ACETYLENE, OXYGEN  FLAMES,  2.7  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results:  O-  CjHj,  Ref.  30,  31;  •  -  Soot,  Ref.  30,  31;  *  -  Soot,  Ref.  32; 
□  -  Temperature  from  AeroChem.  Solid  lines,  no  symbols  -  Ions,  Ref.  33. 
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FIGURE  13  EQUILIBRIUM  BASED  ON  C600H60  AS  SOOT 
TOLUENE/AIR  FLAMES,  101  kPa,  as  in  Figure  11 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results  (Ref.  20):  *  -  Soot,  N2-rich  air;  •  -  Soot,  air;  ▼  -  Soot,  02-rich  air. 

Compare  Figure  11. 
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